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Abstract 

Remarkable success gained by various thermal and statistical approaches in describing the particle 
multiplicities and their ratios has emphasized the formation of a fireball consisting of chemically 
equilibrated hot and dense hadron gas (HG) produced in the ultrarelativistic heavy-ion collisions. In 
an earlier paper referred as I, we proposed a thermodynamically consistent excluded-volume model 
for the HG fireball and we noticed that the model gives a suitable description for various properties 



of multiparticle production and their ratios in the entire range of temperatures and baryon densities. 



Furthermore, a numerical calculation indicates that the model respects causality and the values of 
the transport coefficients (such as shear viscosity to entropy (n/s) ratio, and the speed of sound 
etc.) suitably match with the predictions of other HG models. The aim in this paper is to obtain 
the variations of freeze-out volume in a slice of unit rapidity i.e. dV/dy as well as total volume of 
the fireball with respect to center-of-mass energy (V-SWiv) and confront our model calculations with 
the corresponding thermal freeze-out volume obtained from the Hanbury-Brown-Twiss (HBT) pion 
interferometry. We also test the validity of our model in extracting the total multiplicities as well 
as the central rapidity densities of various hadrons such as 7r + , K + , p, (ft, A, A etc. and 

in getting the rapidity as well as transverse momentum distributions of various particles produced 
in different heavy-ion collider experiments in order to examine the role of any hydrodynamic flow 
on them by matching our predictions with the experimental results. Furthermore, we extend our 
analysis for the production of light nuclei, hypernuclei and their antinuclei over a broad energy 
range from Alternating Gradient Synchrotron (AGS) to Large Hadron Collider (LHC) energies, for 
which the experimental data have started appearing. 
PACS numbers: 12.38.Mh, 12.38.Gc, 25.75.Nq, 24.10.Pa 
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I. INTRODUCTION 



The ultimate goal of ultra-relativistic heavy-ion collisions is to produce a highly excited 
and dense matter which possibly involves a phase transition from a hot, dense hadron gas 
(HG) to a deconfined quark matter called as quark-gluon plasma (QGP) [1-3]. In this state 
the degrees of freedom are those of quarks and gluons only. One hopes that by colliding 
heavy nuclei one can create a fireball with an extremely large energy density extending over 
a sufficiently large space-time volume so that an equilibrated quark-gluon plasma may be 
formed. However, experimental and theoretical investigations made so far reveal that it is 
indeed difficult to get an unambiguous evidence for QGP formation. It is very important 
to understand the dynamics of the collisions in order to suggest any unique signal for QGP. 
Such information can be obtained by analyzing the properties of various particles which are 
emitted from various stages of the collisions. The global properties and dynamics of later 
stages can be best studied via hadronic observables such as hadron yields, ratios, rapidity 
distributions and transverse mass spectra [4]. 

In our previous paper henceforth referred as I [5] , we have proposed a thermodynamically 
consistent excluded-volume model and studied different thermodynamical properties of HG 
such as number density, energy density, pressure etc.. It is indeed surprising that the predic- 
tions of our geometrical model regarding the detailed features of various hadron ratios and 
their variations with respect to the center-of-mass energy (y / ijv r /v) are successfully tested 
with the available experimental data upto Relativistic Heavy-Ion Collider (RHIC) energies. 
We know that if a system passes through a mixed phase in a first order phase transition, 
then it can possess a much larger space-time extent than what we expect from a system if it 
remains in the hadronic phase only. Thus the freeze-out volume V of the fireball formed in 
heavy-ion collisions, is a significant quantity which can give a hint regarding the occurrence 
of mixed phase during a phase transition. Furthermore, the freeze-out volume throws light 
on the collective expansion and/or anisotropic flow of the fluid existing in the fireball before 
a final thermal freeze-out occurs. In this paper, we calculate the variations of chemical 
freeze-out volumes for 7r + , K + , and K~ with respect to y/SNN an d compare our predictions 
with the HBT interferometry data for the thermal freeze-out volume. It gives an indepen- 
dent confirmation pertinent to correctness of the geometrical assumptions involved in our 
model on the space dynamics of the fireball. We also attempt to show the variations of total 



multiplicities as well as the multiplicity in a slice of unit rapidity at midrapidity (y = 0) for 
various emitted hadrons, e.g., ir + , K + , K~, 0, A, S~,fi~, and A. Further, we extend our 
model to deduce the rapidity as well as transverse mass spectra of hadrons and compare 
them with the experimental data available in order to illustrate the role of hydrodynamic 
flow of the fluid. 

In heavy-ion collisions, rapidity densities of produced particles are strongly related to 
the energy density and/or entropy density created in the collisions [6]. Further, the depen- 
dence of transverse mass spectra of hadrons on y / J/vjv can yield insight into the evolution 
of a radial flow present in the dense fluid formed in the collision [7] . Earlier different types 
of approaches have been taken into account for studying rapidity distributions and trans- 
verse mass spectra [8-13]. Hadronic spectra from purely thermal models usually reveal an 
isotropic distribution of particles [8] and hence the rapidity spectra obtained with the ther- 
mal models do not reproduce the features of the experimental data satisfactorily. Similarly, 
the transverse mass spectra from the thermal models reveal a more steeper curve than that 
observed experimentally. The comparisons illustrate that the fireball formed in heavy-ion 
collisions does not expand isotropically in nature and there is a prominent input of col- 
lective flow in the longitudinal and transverse directions which finally causes anisotropy in 
the rapidity and transverse mass distributions of the hadrons after the freeze-out. A lot of 
work has been done in recent years on thermal model calculations incorporating the effect 
of flow [14-23]. Here we mention some kinds of models of thermal and collective flow used 
in the literature. Hydro dynamical properties of the expanding fireball have been initially 
discussed by Landau and Bjorken for the stopping and central-rapidity regimes, respectively 
[9]. However, collisions even at RHIC energies reveal that they are neither fully stopped, 
nor fully transparent. As the collision energy increases, the longitudinal flow grows stronger 
and leads to a cylindrical geometry as postulated in Ref. [10-12]. They assume that the 
fireballs are distributed uniformally in the longitudinal direction and demonstrate that the 
available data can consistently be described in a thermal model with inputs of chemical 
equilibrium and flow, although they have used the experimental data for small systems only. 
They use two simple parameters : transverse flow velocity (j3 r ) and temperature (T) in their 
models. In Ref. [13], non-uniform flow model is used to analyze the spectra specially to 
reproduce the dip at midrapidity in the rapidity spectra of baryons by assuming that the 
fireballs are distributed non-uniformly in the longitudinal phase space. In Ref. [24-26], the 



rapidity-dependent baryon chemical potential has been invoked to study the rapidity spectra 
of hadrons. In certain hydrodynamical models [27,28], measured transverse momentum (pr) 
distributions in Au — Au collisions at y/s NN = 130 GeV [29-31] have been described suc- 
cessfully by incorporating a radial flow. In Ref. [32], rapidity spectra of mesons have been 
studied using viscous relativistic hydrodynamics in a 1+1 dimension assuming a non-boost 
invariant Bjorken flow in the longitudinal direction. They have also analyzed the effect of 
the shear viscosity on the longitudinal expansion of the matter. Shear viscosity counteracts 
the gradients of the velocity field, as a consequence slows down the longitudinal expansion. 
In this paper, we attempt to calculate the rapidity density of various particles at midra- 
pidity by using our excluded-volume model of HG and we find a good agreement between 
our model calculations and the experimental data. However, when we calculate the rapidity 
distributions of various particles, we find that the distribution always takes a narrow shape 
in comparison to the experimental data and thus indicating that our thermal model alone is 
not capable to describe the experimental data in forward and backward rapidity regions. In 
a similar way, the transverse mass spectra of hadrons as obtained in our model again does 
not describe the experimental data properly. Our analysis clearly necessitates the presence 
of an additional flow factor to be introduced in the model in order to provide a suitable 
description of the data. 

We plan to study the rapidity distributions and transverse mass spectra of hadrons using 
our model in I [5] for a hot, dense HG which is based on an excluded-volume approach [33] 
formulated in a thermodynamically consistent way. We have assigned an equal hard-core 
size to each type of baryons in the HG while the mesons are treated as pointlike particles. 
We use the chemical freeze-out criteria to relate the thermal parameters T and /ie with the 
collision energies [5] . We further emphasize that we have used this hadronic equation of state 
(EOS) together with a suitable EOS for QGP in order to obtain the QCD phase diagram and 
the associated critical point [34]. The plan of the paper runs as follows : the next section 
is dedicated to the formulation of our model for the study of the rapidity distributions 
and transverse mass spectra of hadrons using purely thermal source. In the third section, 
we modify the formula for the rapidity distributions by incorporating a flow velocity in 
the longitudinal direction and similarly the formula for the transverse mass spectra is also 
modified by incorporating a collective flow in the longitudinal as well as in the transverse 
direction. In section IV, we compare the experimental data with our predictions regarding 



rapidity density, transverse mass spectra and transverse momentum spectra of hadrons at 
RHIC energies (130 GeV and 200 GeV) and LHC energy (2.76 TeV). We have also shown 
the rapidity spectra of 7r~ at Super Proton Synchrotron (SPS) energies. As v /saT/v increases, 
the longitudinal as well as the transverse flow also increase which means that the experiments 
will observe a larger amount of collective flow at LHC. We have also deduced total mean 
multiplicities, mid-rapidity yields of various hadrons and their rapidity densities and studied 
their variations with v /J/vjv. Thus we determine the freeze-out volume and/or dV/dy in order 
to ascertain whether the emissions of all hadrons occur from the same hypersurface of the 
fireball. We also analyze the experimental data on the production of light nuclei, hypernuclei 
and their antiparticles over a broad energy going from AGS to RHIC energies. This analysis 
thus fully demonstrates the validity of our EOS in describing all the features of a hot, dense 
HG. Finally, in the last section we succinctly give the conclusions and summary. 



II. HADRONIC SPECTRA WITH THE THERMAL SOURCE 



A. Rapidity Distributions 

To study the rapidity distributions and transverse mass spectra of the hadrons, we extend 
excluded- volume model as proposed in I [5], which was found to describe the hadron ratios 
and yields at various v /s/v r /v from lower AGS energies upto RHIC energies with remarkable 
success. We have incorporated the excluded-volume correction directly in the grand canon- 
ical partition function of HG in a thermodynamically consistent manner. We obtain the 
number density nf for ith species of baryons after excluded-volume correction as follows [5]: 

nT = (1 - R)h\ - h\ 2 ^ + A, 2 (l - R)li (1) 
where R = n^V® is the fractional occupied volume by the baryons [5,34]. V t ° = Air r 3 /3 

i 

is the eigen-volume of each baryon having a hard-core radius r and A« is the fugacity of ith 
baryon. Further, Jj is the integral of the baryon distribution function over the momentum 
space [5]. 

We can rewrite Eq. (1) in the following manner : 
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This reveals that the invariant distributions are [10,11] : 



d 3 N t 9i V\ 



dk 3 (2tt) ; 



-Ai(l - i?) 



( r / ; ( y ) A, 



Ei 



(3) 



/•; ( \ 2 
earp ( y J + A* 



Using : 



E, 



d 3 N t 



dNi 



dk 3 dy rriT drriT d(j> ' 



(4) 



we get : 



dNi 



= 9iV\ 

dy rriT drriT d(p (2n) 3 



l-R)- A 



dR \ 



Ei 



w 

-Ai(l -R)- 



exp 



T 



+ \i 



Ei 



exp 



E\ I 2 

yl+A, 



(5) 



Here y is the rapidity variable and transverse mass tut = \/m 2 + p^ 2 . Also Ei is the energy 
of ith baryon and V is the total volume of the fireball formed at chemical freeze-out and iVj 
is the total number of ith baryons. We assume that the freeze-out volume of the fireball for 
all types of hadrons at the time of the homogeneous emissions of hadrons remains the same. 

By inserting E^ = m T coshy in Eq. (5) and integrating the whole expression over trans- 
verse component we can get the rapidity distributions of baryons as follows: 
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Eq.(6) gives the rapidity distributions of baryons arising due to a stationary thermal 
source. Here we evaluate the rapidity distributions of hadrons at px = as the data have 
been taken for the p T = case. 

Similarly, the rapidity density of mesons can be had by using the following formula : 
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Here g m , \ m are the degeneracy factor and fugacity of the meson m, respectively. 



B. Transverse Mass Spectra 

We use Boltzmann statistics in deriving formula for the transverse mass spectra because 
we want to calculate spectra of hadrons only at RHIC energies, where the effect of quantum 
statistics is found to be negligible [10]. In the Boltzmann's approximation, Eq.(5) can be 
reduced to a simple form : 
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Putting Ei = m T coshy in Eq.(8) and integrating over rapidity (y) we get the transverse 
mass spectra as follows : 
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where K x ( -=- ) is the modified Bessel's function and is given by the following expression 



Similarly transverse mass spectrum of mesons can be evaluated as follows : 
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III. HADRONIC SPECTRA WITH THE EFFECT OF FLOW 

In the previous section, we have obtained the expression for rapidity distributions as well 
as transverse mass spectra arising from a stationary thermal source only. In this section, we 
modify the expression for rapidity spectra i.e. Eq. (6), by incorporating a flow velocity to a 
stationary thermal system in the longitudinal direction and thus we attempt to explain the 
experimental data in the full rapidity region. The resulting rapidity spectrum of ith hadron, 
after the incorporation of the flow velocity in the longitudinal direction is [10,11]: 
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where ( 1 -\ can be calculated by using Eq.(6) for the baryons and by Eq.(7) for the 

\ dy J th 

mesons. The average longitudinal velocity is [13,35]: 



{(5 L )=tanh{^f). (14) 

Here i] max . is an important parameter which provides the upper rapidity limit for the longitu- 
dinal flow velocity at particular ^s^jy and it's value is determined by the best experimental 
fit. The value of r] max , increases with the increasing ^s^n and hence (3l also increases. 

In the case of transverse mass spectra, we incorporate flow velocity in both the directions, 
longitudinal as well as transverse directions in order to describe the experimental data 
satisfactorily. However, we assume radial type of flow velocity in the transverse direction 



which imparts a radial velocity boost on top of the thermal distribution. We define the four 
velocity field in both the directions as follows [36] : 



u^(p,n) = (coshp coshr),e r sinhp,coshp sinhrj), 
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where p is the radial flow velocity in the transverse direction and rj is the flow velocity in 
the longitudinal direction. Once we have defined the flow velocity field, we can calculate the 
invariant momentum spectrum by using the following formula [10,37] : 
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In the derivation of Eq.(16) we assume that an isotropic thermal distribution of hadrons is 
boosted by the local fluid velocity w M . Now the resulting spectrum can be written as follows 
[10]: 
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The freeze-out hypersurface d<7\ in Eq.(16) is parametrized in cylindrical coordinates (r, 0, (), 
where the radius r can lie between and R i.e. the radius of the fireball at freeze-out, the 
azimuthal angle lies between and 2n, and the longitudinal space-time rapidity variable 
( varies between —r] max and r\ max . Now, integrating Eq.(17) over <p as well as (, we get the 
final expression for the transverse mass spectra [10] : 
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is the modified Bessel's function given by : 



Pt sinhp\ 



T 



-) 



2w 

27 1 exr 



Pt sinhp cos<p 
T 



(19) 



where p is given by p = tanh 1 f3 r , with the velocity profile chosen as (3 r = (3 S [10,11]. 
(3 S is the maximum surface velocity and is treated as a free parameter and £ = (r/R^j. The 
average of the transverse velocity can be evaluated as [31]: 




FIG. 1: Energy dependence of the freeze-out volume for the central nucleus-nucleus collisions. 
The symbols are the HBT data for freeze-out volume Vhbt for the n + [38]. A' ,B' and C are the 
total freeze-out volume and A ,B and C depict the dV/dy as found in our model for ir + , K + and 
K~ , respectively. 



In our calculation we use a linear velocity profile, n — 1 and R is the maximum radius 
of the expanding source at freeze-out (0 < £ < 1) [31]. Similarly following equation can be 
used to calculate transverse mass spectra for mesons : 
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IV. RESULTS AND DISCUSSIONS 

We have used freeze-out temperature and baryon chemical potential as determined by 
fitting the particle ratios as described in I [5]. Then we use suitable parametrizations for T 
and ix b in terms of center-of-mass energies. We have used hard-core radius r = 0.8 fm for 
all types of baryons. We have considered all the particles and the resonances upto mass of 




FIG. 2: Variations of total multiplicities of tt + , K + , K~ , <fi, A, S _ , (Q~ + and A with 

respect to center-of-mass energy predicted by our model. Experimental data measured in central 
Au — Au/Pb — Pb collisions [39-55] have also been shown for comparison. In this figure, A, B, 
C, D, E, F, G, and H represent the multiplicities of ir + , K + , K~, 4>, A, E~, + and A, 
respectively. 

2 GeV/c 2 in our calculation. We have used resonances having well denned masses, widths 
etc. and branching ratios for sequential decays are also suitably incorporated. 

Although the emission of hadrons from a statistical thermal model essentially invokes the 
idea of equilibration, it does not reveal any information regarding the existence of a QGP 
phase before hadronization. However, if the constituents of the fireball have gone through 
a mixed phase, volume V of the fireball at freeze-out is expected to be much larger than 
what we expect from a system if it remains throughout in the hadronic phase only. In 
Fig. 1, we have shown V and dV/dy as obtained in our excluded- volume model and their 
variations with the center-of-mass energy. For this purpose, we have used the data for total 
multiplicities of 7r + , K + , K~ and after dividing with their corresponding number densities 
obtained in our model, we get V for n + , K + , and K~, respectively. Similarly for deducing 
dV/dy, we use the data for dN/dy and divide by the corresponding number density. We 
have compared predictions of our model with the experimental data obtained from the pion 
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FIG. 3: Variation of rapidity distributions of various hadrons with respect to y / ^s/v r iv at midrapidity. 
Lines show our model calculation. Symbols are the experimental data [7,56,57]. 
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FIG. 4: Energy dependence of rapidity distributions of hadronic species at midrapidity (y = 
0) for central nucleus-nucleus collisions. Lines are the model calculations and symbols are the 
experimental results [7,56,57] at various ^snn- 



interferometry (HBT) [38] which in fact reveal thermal (kinetic) freeze-out volume. Our 
results support the finding that the decoupling of strange mesons from the fireball takes 
place earlier than the 7r-mesons. Moreover, a flat minimum occurs in the curves around 
energy 8 GeV, then the volume rapidly increases. At RHIC energy, the volume of the 
fireball is around 12000 fm? which is much larger than the volume of gold nucleus used in 
the collisions. 

In order to calculate total multiplicities of hadrons, we first determine the total freeze-out 
volume for K + by dividing the experimentally measured multiplicities of K + at different 
center-of-mass energies with it's number density as calculated in our model at different 
center-of-mass energies. We assume that the fireball after expansion, achieves the stage of 
chemical equilibrium and the freeze-out volume of the fireball for all types of hadrons at the 
time of the homogeneous emissions of hadrons remains same for all particles. This freeze-out 
volume thus extracted for K + , has been used to calculate multiplicities of all other hadrons 
from corresponding number densities at different y / s/v r /v. Figure 2 shows the center-of-mass 
energy dependence of multiplicities of hadrons n + , K + , K~, 0, A, S~, (fi - + and 
A as predicted by our model calculation. We also show here corresponding experimental 
data measured in central Au — Au/Pb — Pb collisions [39-55] for comparison. We observe 
an excellent agreement between our model calculation and experimental data for the total 
multiplicities of all particles except 0, S~, and Q~ , where we see small deviations. Again, 
the thermal multiplicities for all these particles are larger than the experimental values. This 
analysis again suggests a new and different mechanisms for the production of these particles. 

In Fig.3, we compare the experimental midrapidity data [7,56,57] of various hadron species 
over a broad energy range from AGS to RHIC energies with the results of our model calcu- 
lations. We use the same freeze-out volume of the fireball for each hadron as extracted for 
K + in our model calculation. We also show the comparison between the results calculated 
using our model with the Boltzmann's statistics, as used in our previous paper [58], with 
the results calculated with full quantum statistics. As expected, both results differ only at 
lower energies. However, the results with quantum statistics give much closer agreement 
with the experimental data. Figure 4 shows the energy dependence of midrapidity distri- 
butions of various hadrons like p, 0, p, S~, and S+. We observe that the results obtained 
in our model show a close agreement with the experimental data [7,56,57]. However, our 
model calculation again differs at higher energies for multistrange particles in comparison 
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FIG. 5: Rapidity distribution of ir + at ^/snn = 200GeV. Dotted line shows the rapidity distribu- 
tion calculated in our thermal model. Solid line and dashed line show the results obtained after 
incorporating longitudinal flow in our thermal model. Symbols are the experimental data [59]. 

to experimental data, which appears to be a common problem of all thermal models. 

In Fig.5, we present the rapidity distribution of 7r + for central Au + Au collisions at 
y/sNN = 200 GeV over full rapidity range. Dotted line shows the distribution of ir + due to 
stationary thermal source which describes only the midrapidity data while it fails to describe 
the experimental data at other rapidities. Solid line shows the rapidity distributions of n + 
after incorporation of longitudinal flow in our thermal model and again in a good agreement 
with the experimental data [59]. After fitting the experimental data, we get the value of 
Vmax. = 3.20 and hence the longitudinal flow velocity (5l = 0.922 at y/s NN = 200 GeV. 
For comparison and testing the appropriateness of parameter, we also show the rapidity 
distributions at different values of the parameter, r) max , = 2.80 i.e. (3l = 0.88, which is 
shown by the dashed line in the figure. We find that the results now differ. 

In Fig.6, we analyze the experimental data [60] on rapidity distributions of 7r~ at various 
t/snn in terms of our model calculation to constrain the allowed distribution of the longitu- 
dinal velocities of the fluid elements. We find that our model results are in good agreement 
with the experimental data at these energies. Figure 7 demonstrates the variation of longi- 
tudinal flow velocities extracted in our model with respect to y^S/vjv and it shows that the 
longitudinal flow velocity increases with the increasing ^snn- We compare our results with 




FIG. 6: Rapidity distribution of ir at y / ijv r /v = 8.76 GeV, 12.3 GeV and 17.3 GeV. Lines are our 
model calculation and symbols are experimental data [60]. 
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FIG. 7: Variations of the longitudinal flow velocities with respect to ^snn- Open circles are the 
results of Ref . [35] . 

that calculated in Ref. [35] and find that a slight difference exists between these two results. 

Figures 8 and 9 demonstrate the rapidity distributions of proton and anti-proton, respec- 
tively at ^snn = 200 GeV. Again, our curves describe successfull the experimental data 
[61]. Similarly in Figures 10 and 11 we show the rapidity distributions for K + and K~ , 
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FIG. 8: Rapidity distribution of proton for central Au — Au collisions at ^Jsnn = 200 GeV . 
Dotted line shows the rapidity distribution due to purely thermal source and solid line shows the 
result after incorporating the longitudinal flow velocity in our thermal model. Symbols are the 
experimental data [61]. 
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FIG. 9: Rapidity distribution of anti-proton for central Au — Au collisions at y / sjv r /v = 200 GeV. 
Dotted line shows the rapidity distribution due to purely thermal source and solid line shows the 
result after incorporating the longitudinal flow in our thermal model. Symbols are the experimental 
data [61]. 
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FIG. 10: Rapidity distribution of K + at ^Jsnn = 200 GeV. Dotted line shows the rapidity distri- 
bution due to purely thermal source and solid line shows the result after incorporating longitudinal 
flow in the thermal model. Symbols are the experimental data [62]. 

respectively at ^s NN = 200 GeV. We have used the same value of parameter j3 L and V for 
each particle. 

Figure 12 shows the analysis of the transverse mass spectra of ir + and proton for the 
most central collisions of Au — Au at ^snn = 130 GeV in terms of our model calculation. 
We have neglected the contributions from resonance decays in our calculations here because 
the contributions from resonance decays affect the transverse mass spectra only at the lower 
transverse mass side i.e. rriT < 0.3 GeV. So, our model calculation show some difference 
with the experimental data [31] for mr < 0.3 GeV but there is a good agreement between 
our model calculation and experimental results for m T > 0.3 GeV . We also show our 
model results for n + at different values of f3 s just for comparison. Similarly Figure 13 shows 
the transverse mass spectra of K + and for the most central collisions of Au — Au at 
v /sjv r /v = 130 GeV . Our calculation is quite successfull in describing the experimental data 
[31,64]. However, our model fails again in the case of multistrange particle <fi. We need 
to devise some different mechanism to explain the experimental results for multistrange 
particles. By fitting the experimental results, we get the value of /3 S = 0.45 and hence 
transverse flow velocity (3 r = 0.30 at y/s NN = 130 GeV, and is less than that extracted in 
Ref. [31]. It may be due to different freeze-out temperature used in our model. 





FIG. 11: Rapidity distribution of K~ at ^/smn = 200GeV . Dotted line shows the rapidity distri- 
bution due to purely thermal source and solid line shows the result after incorporating longitudinal 
flow in our thermal model. Symbols are the experimental data [63]. 
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FIG. 12: Transverse mass spectra for 7r + and proton for the most central collision at ^/snn = 
130GeV. Dashed and dotted lines are the results of transverse mass spectra, calculated in our 
thermal model, for it + and proton, respectively. Solid and dash-dotted lines are the results for ir + 
and proton, respectively, obtained after incorporation of flow in our thermal model. Symbols are 
the experimental data [31]. 
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FIG. 13: Transverse mass spectra for K + and 4> for the most central collision at ^/snn = 130 GeV. 
Dashed and dotted lines are the results of transverse mass spectra, calculated in our thermal 
model, for K + and <f>, respectively. Solid and dash-dotted lines are the results for ir + and proton, 
respectively, obtained after incorporation of flow in thermal model. Symbols are the experimental 
data [31,64]. 

In Fig. 14, we show the transverse mass spectra of n + and proton for the most central 
collisions of Au + Au at y/SNN = 200 GeV . From this figure it is clear that the spectra 
due to stationary thermal source do not satisfy the experimental data [63] at higher m T 
while the thermal source with collective flow can describe the experimental data very well. 
For comparison, we also show the transverse mass spectra for n + at different j3 s . In Fig. 

15, we present the transverse-mass spectra for K + and at ^snn = 200 GeV. Lines are 
our model calculations and the symbols represent the experimental data [63,64]. Again, our 
model fails to describe the experimental data for at lower m T . At this v /sjvjv the value 
of (5 S) after best fitting, comes out to be 0.50 and hence transverse flow velocity f3 r = 0.33. 
The transverse flow velocity is able to reproduce the transverse mass spectra of almost all 
the hadrons at y / sJvjv = 200 GeV. We also see that the transverse flow velocity increases 
with the increasing ^snn- 

We also present the pr- spectra of hadrons at various center-of-mass energies. In Fig. 

16, we show the (p T ) spectra of 7T + , K + and p for the most central collisions of Au — Au 
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FIG. 14: Transverse mass spectra for tt + and proton for the most central collisions at y/sJm = 
200 GeV. Dashed and dotted lines are the transverse mass spectra due to purely thermal source 
for 7r + and proton, respectively. Solid and dash-dotted lines are the results for ir + and proton, 
respectively obtained after incorporation of flow in thermal model. Symbols are the experimental 
data [63]. 

at y/SNN = 200 GeV. Our model calculations are in close agreement with the experimental 
data [63]. In Fig. 17, we show the p T spectra of 7r~, K~ and p for the Pb — Pb collisions at 
y/SNN — 2.76 TeV at LHC. Our model calculations give a good fit of the experimental results 
[65]. We also compare our results for p spectrum with the hydrodynamical model [66], which 
explains successfully the n~ , and K~ spectra but strongly fails to describe the p spectrum. 
We find that our model results are closer with the experimental data than that of Ref.[66]. 
In Ref. [66], Shen et al. have employed 2 + 1-dimensional viscous hydrodynamics with the 
lattice QCD-based EOS. They use Cooper-Frye prescription to implement kinetic freeze- 
out by converting the hydrodynamic output to particle spectra. Due to lack of a proper 
theoretical and phenomenological knowledge, they use the same parameters for Pb — Pb 
collisions at LHC energy, which was used for Au — Au collisions at v /saT/v = 200 GeV . 
Furthermore, they use the temperature independent i]/s ratio in their calculation. After 
fitting the experimental data, we get (5 S = 0.80 and hence f3 r = 0.53 at this energy which 
indicates that the collective flow is stronger at this energy than that observed at RHIC 
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FIG. 15: Transverse mass spectra for K + and cf> for the most central collision at ^Jsnn = 200GeV. 
Dashed and dotted lines are the transverse mass spectra due to purely thermal source for K + and 
4>, respectively. Solid and dash-dotted lines are the results for K + and (f>, respectively, obtained 
after incorporation of flow in thermal model. Symbols are the experimental data [63,64]. 
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FIG. 16: Transverse momentum spectra for ir + , p, and K + for the most central Au — Au collision at 
^Jsnn = 200 Gel/. Lines are the results of our model calculation and symbols are the experimental 
results [63]. 
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FIG. 17: Transverse momentum spectra of various hadrons for the most central collisions of Pb—Pb 
at ^/iivjv = 2.76 TeV from LHC. Lines are the results of model calculations and symbols are the 
experimental results [65]. Thick-dashed line is the prediction of viscous- hydrodynamical model [66] 
for p. Dashed line is the prediction of our model calculation for cj) meson. 




FIG. 18: Energy dependence of the ifi/ir, (n = 1.5(7r + +vr - )) ratio at midrapidity. Dash-dotted line 
is the result of statistical hadron gas model (HGM) [67]. Solid line represents our model calculation 
and dashed line is the result of UrQMD model [68] . Symbols are the experimental data points [7] . 




FIG. 19: The energy dependence of anti-baryon to baryon yield and antinuclides to nuclides ratios. 
Lines are our model calculation and symbols represent the experimental data [74]. 

energies. We also predict the pt spectra for <fi meson at this energy. For comparison, we 
also show the spectra at different values of parameter j3 s = 0.88. 

In order to emphasize how our thermal model fails for multistrange particles, we show in 
Fig. 18 the variations of the <j>/n, (it = 1.5(7r + + ratio at midrapidity with -y/sjvjv. We 
compare our results with the experimental data [7] and statistical hadron gas model (HGM) 
[67]. We find that our model provides a good description of the experimental data at lower 
energies but fails to describe the data at higher energies. HGM model, which has included 
the strangeness saturation factor 7 S to account for the partial equilibration of strangeness, 
describes the data well at the lower and upper energies while it completely disagrees with 
the experimental data at the intermediate energies. We also show the results of the hadronic 
transport model Ultra-relativistic Quantum Molecular Dynamical (UrQMD) [68] in which 
kaon coalescence mechanism is used for the production of 0. The predictions of UrQMD lie 
well below the experimental data at all energies which suggests that this mechanism fails 
badly in explaining the production mechanism of 0. 

We present an analysis of light nuclei, hypernuclei and their anti-particles using the 
chemical freeze-out concept within our model calculation. The productions of light nuclei 
and hypernuclei at chemical freeze-out points may not be appropriately calculated, because 
their binding energies are of the order of few MeV and the chemical freeze-out temperatures 




FIG. 20: The energy dependence of various baryons, antibaryons, nuclei, and antinuclei yield ratios. 
Lines are our model calculation and symbols represent the experimental results [74]. 

are around 100 — 165 MeV. But we know that the relative yield of particles composed of 
nucleons is mainly determined by the entropy per baryon which is fixed at chemical freeze- 
out line in our model [5]. This was first outlined in [69] and was subsequently emphasized 
in [70]. This constitutes the basis of thermal analyses for the yields of light nuclei [71,72]. 
Thus the production yields of light nuclei and hypernuclei is entirely governed by the entropy 
conservation. Recently the first measurement of the lightest (anti) hypernuclei was done by 
the STAR experiment at RHIC [73]. We want to ask an interesting question whether the 
productions of antinuclei can also be explained by our model. In this paper, we predict and 
throw light on the production yields of light nuclei, hypernuclei, and heavy baryons (anti- 
baryons) within our thermal model approach and we compare our thermal model calculations 
with the experimental data. The thermal freeze-out parameters for these yields are also taken 
as the same as used for other ratios. 

In Fig. 19, we show the energy dependence of p/p, A/A, He 3 /He 3 , He 4 / He 4 and d/d 
yield ratios over a very broad energy range. We take three times of volume and mass of each 
nucleon to calculate the volume and mass of He 3 in our calculation. Similarly, we calculate 
densities of other nuclei and hypernuclei using the similar method. We compare our results 
with the experimental data [74] and find a good agreement between these two. We stress that 
the agreement between the experimental and the calculated curves for the ratio He 3 /He 3 



is a powerful argument that indeed entropy conservation governs the production of light 
nuclei. In Fig. 20, we show the variations of A/p,d/p with ^snn- We compare our results 
with the experimental data points [74]. We also predict the energy dependence of nuclei as 
well as antinuclei yield ratio such as H\/He 3 and H^/He 3 over a broad energy range going 
from lower AGS energies to LHC energy. Some data have started appearing and we find 
that our model gives a good fit to the data. 

V. CONCLUSIONS 

In conclusion, we find that our model provides a good fit to the variations of total multi- 
plicities as well as mid-rapidity densities of various particles and we deduce a large freeze-out 
volume of the fireball at RHIC energy and this picture supports the idea of a mixed phase 
after QGP formation before hadronization because a huge size of a homogeneous fireball 
source can only arise if a mixed phase has occurred before the formation of a hot, dense 
HG. Further, we present an analysis of rapidity distributions and transverse mass spectra 
of hadrons in central nucleus-nucleus collisions at various center-of-mass energies using our 
recently proposed equation of state (EOS). We see that the stationary thermal source alone 
cannot describe the experimental data fully unless we incorporate flow velocities in the lon- 
gitudinal as well as in the transverse direction in our thermal model and as a result our 
modified model predictions show a good agreement with the experimental data. Our anal- 
ysis shows that a collective flow develops at each ^s^n which increases further with the 
increasing tJsnn- The description of the rapidity distributions and transverse mass spectra 
of hadrons at each ^snn matches very well with the experimental data. Although, we are 
not able to describe successfully the spectra for multistrange particles which suggests that 
a somewhat different type of mechanism is required. We find that the particle yields and 
ratios measured in heavy-ion collisions are described well by our thermal model and show an 
overwhelming evidence for chemical equilibrium at all beam energies. The rapidity distribu- 
tions and transverse mass spectra which essentially are dependent on thermal parameters, 
also show a systematic behaviour and their interpretations most clearly involve the pres- 
ence of a collective flow in the description of the thermal model. We have also found that 
our model together with a quasiparticle EOS for QGP gives the QCD phase boundary and 
the location of critical point almost overlapping on the freeze-out curve as predicted earlier 



[34]. In conclusion, we find that our model provides an excellent description for the thermal 
and hadronic properties of the HG fireball and associated phase transitions existing in the 
ultra-relativistic heavy-ion collisions. 

In summary, we have formulated a thermodynamically consistent EOS for a hot and 
dense HG fireball created in the ultra-relativistic heavy-ion collisions. We have incorporated 
a geometrical hard-core size for baryons and antibaryons only. In our description mesons may 
possess a size but they can penetrate and overlap on each other. Our prescription has also 
made use of full quantum statistics [5]. It is encouraging to note that our results for particle- 
ratios and their energy-dependence, total multiplicities, rapidity densities and rapidity and 
momentum spectra of various particles consistently match with the experimental data. We 
should emphasize that the hadrons emerging from a completely thermalized HG fireball do 
not give any information regarding a primordial QGP phase existing before HG. However, 
the large freeze-out volume of the fireball obtained in this analysis hints at a mixed phase 
occurring before pure HG phase. 
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